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ABSTRACT 

The globular cluster NGC 6388 is among the most massive clusters in our Milky Way and has been the subject of many studies. 
Recently, two independent groups found very different results when measuring its central velocity-dispersion profile with different 
methods. While Lutzgendorf et al. (2011) found a rising profile and a high central velocity dispersion (23.3 km s“*), measurements 
obtained by Lanzoni et al. (2013) showed a value 40% lower. The value of the central velocity dispersion has a serious impact 
on the mass and possible presence of an intermediate-mass black hole at the center of NGC 6388. The goal of this paper is to 
quantify the biases arising from measuring velocity dispersions from individual extracted stellar velocities versus the line broadening 
measurements of the integrated light using new tools to simulate realistic IFU observations. We use a photometric catalog of NGC 
6388 to extract the positions and magnitudes from the brightest stars in the central three arcseconds of NGC 6388, and create a 
simulated SINEONI and ARGUS dataset. The construction of the IFU data cube is done with different observing conditions (i.e., Strehl 
ratios and seeing) reproducing the conditions reported for the original observations as closely as possible. In addition, we produce an 
N-body realization of a ~ 10® Mq stellar cluster with the same photometric properties as NGC 6388 to account for unresolved stars. We 
find that the individual radial velocities, i.e. the measurements from the simulated SINEONI data, are systematically biased towards 
lower velocity dispersions. The reason is that due to the wings in the point spread function (PSF) of adaptive optics (AO) corrected 
data sets, the velocities get biased towards the mean cluster velocity. This study shows that even with AO supported observations, 
individual radial velocities in crowded fields do not reproduce the true velocity distribution. The ARGUS observations do not show this 
kind of bias but were found to have larger uncertainties than previously obtained. We find a bias towards higher velocity dispersions in 
the ARGUS pointing when fixing the extreme velocities of the three brightest stars but find those variations are within the determined 
uncertainties. We rerun Jeans models and fit them to the kinematic profile with the new uncertainties. This yields a black-hole mass 
of M. = (2.8 ± 0.4) X 10'' Mq and M/L ratio M/L = (1.6 ± Q.\)MqIL i^, consistent with our previous results. 
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1. Introduction 

Measuring central velocity-dispersion profiles in Galactic glob¬ 
ular clusters aims to detect the crucial rise that would reveal the 
presence of an intermediate-mass black hole (IMBH) in the cen¬ 
ter. In the past years, several techniques have been established to 
derive the central velocity-dispersion pr ofile in globular cluste rs. 
Starting with long-slit observations (e.g. lPeterson et aOl989l) to 
measure individual velocities and integrated light, integral field 
units (IFUs ) soon became a valu able tool to measure the central 
kinematics. iDubath et al.l (Il994t) performed extensive numerical 
simulations of the velocity-dispersion determination from inte¬ 
grated light for the globular cluster M15. They found that inte¬ 
grated light measurements over small areas (~ 1 arsec“) suffer 
from large statistical errors due to shot noise of a few massive 
stars, but that a larger coverage of the region of interest repro¬ 
duces well the underlying velocity distribution. Besides using 
the integrated spectra of an IFU for the velocity dispersion deter¬ 
mination, there is also the possibility of obtaining individual ve¬ 
locities of the stars in the field-of-view. This has been done using 


adaptive optics (AO) supported observations and extracting the 
spectr a from the central spaxel of each star (e.g. iLanzoni et al.l 
12013l) or by extracting individual spe ctra using a sophistic ated 
deconvolution technique developed bv iKamann et alJ(|2013l) . 

The first detection of an IMBH in a globular clu st er us- 
ing IFU spectroscopy was reported by iNovola et alJ (l2008l) 
who studied the central region of the globular cluster u) 
Centauri and found the velocity dispersion to be rising to- 
wards the center. This det e ction, however, was challenged by 
Ivan der Marel & AndersonI (120101) who did not see such a rise 
in their velocity-dispersion profile when using proper motions 
and a different photometric center. Aft er this, many IFU obser¬ 
vations of globular clusters followed (iLiitzgendorf et al.]l201 iL 
l20 1 31 : l^ldmeier et al Jl20 1 31) which showed signatures of IMBH 
in the center in some Galactic globular clusters. However, con- 
tradic tory results using in dividual radial vel ocities instead of 
IFUs (ILanzoni et al.ll20l'^ or proper motions (iMcNamara et al.l 
12012l) kept the discussions going on whether or not IMBHs in 
the centres of globular clusters exist. Furthermore, X-ray and 



































2 


Nora Liitzgendorf et al.: A re-evaluation of the central velocity-dispersion profile in NGC 6388 


radio observations of the central regions of globular clusters re¬ 
sult in no detection and low upper limits on the black-hole mass 
(e.g. lMaccarone et al.l2005HStrader et al.l2012h . To date, there is 
no evidence of accreting IMBHs in globular clusters. The upper 
limits on black-hole masses from these observations are, how¬ 
ever, very dependent on the knowledge of gas densities and ac¬ 
cretion variability. 

The case of NGC 6388 is an excellent ex ample of contradict¬ 
ing me asurements with different techniques. iNovola & Gebhardtl 
(l2006ll found a shallow cusp in the central region of the surface 
brightness profile of NGC 6388. N-body simulations showed 
that this is expected for a cluster hosting an intermediate-mas s 
black hol e (iBaumgardt et al.ll2005 jNovola & Baumgarddl20ini . 
However, IVesperini & Trentil (1201 Oh argued that shallow cusps 
can also be produced by clusters that are in the process of core 
collapse and do not exhibit IMBHs, and the interpretation of 
the cusp is therefore still under debate. iLanzoni et al.l (l2007h in¬ 
vestigated the projected density profile and the central surface 
brightness profile with a combination of HST high-resolution 
and ground-based wide-held observations. They found the ob¬ 
served prohles are w ell reproduced by a multimass, isotropic, 
spherical iKingl d 19661) model, with an added central black hole 
with a mass of ~ 5.7 x IO^Mq. Our group pres ented the hrst cen¬ 
tral kin ematic measurements of NGC 6388 in lLiitzgendorf et ^ 
(1201 ih and found the velocity dispersion prohle rising in the cen¬ 
ter up to 23.3 km s ' indicating the presence o f an IM BH with 
a mass of ~ 10 "'Mq. However, ILanzoni et alJ (1201 3l) recently 
claimed the velocity dispersion in the center of NGC 6388 to 
be lower by 40%. Using individual radial velocities extracted 
from the adaptive optics supported IFU of the SINFONI instru¬ 
ment, they found the velocity-dispersion to be 13.2 km s * using 
52 individual velocities instead of 23.3 km s * from the ARGUS 
measurements. This high discrepancy is worrisome, as it af¬ 
fects the measurement on the existence and mass of a possible 
intermediate-mass black hole in the center. 

Our group has therefore set out to investigate the effect on 
the velocity dispersion measurement when using individual ve¬ 
locities and integrated light. In Section |2 we introduce the two 
different observations and data sets. Section[3]describes the sim¬ 
ulations that were done to reconstruct the two different IFU ob¬ 
servations and SectionlUpresents the extraction of the kinemat¬ 
ics. The results of the Monte Carlo simulations and a summary 
can be found in Section|5]and|6l respectively. 


2. The data sets 

Two independent data sets provide central velocity dispersions 
for NGC 6388. The first one was observed with the GIRAFFE 
spectrograph of the FLAMES (Eiber Large Array Multi Element 
Spectrograph) instrument at the Very Large Teles cope (VLT) us¬ 
ing th e ARGUS mode (Large Integral Eield Unit. lPasauini et aP 
I 2 OO 2 I) . The second observation covered a smaller area in the 
center using the high spatial resolution cap abilities of SINEONI 
(lEisenhauer et al.lUOojUBonnet et al]|200^ . a near-IR (1.1-2.45 
^m) integral field spectrograph fed by an adaptive optics module 
also mounted on the VLT. In this section we shortly describe the 
two observations and their conditions. 

The ARGUS observations were performed during two nights 
(2009-06-14/15, ESO proposal ID: 083.D-0444; PI: Noyola) 
with an average seeing of 0.8" (EWHM). The lEU unit was 
set to the 1 : 1 magnification scale (pixel size: 0.3", 14 x 22 
pixel array) with the LR8 filter (820 - 940 nm, R = 10400) and 
pointed to three different positions to cover the central area of 
NGC 6388. The kinematics were obtained from the analysis of 
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Fig-1- Veloci t y disp ersion as a function of magnitude for the 
ILanzoni et al.l (120131) sample. The sample shows a trend of de¬ 
creasing velocity dispersion with increasing magnitude which 
cannot be explained by mass segregation. 



the Calcium Triplet (~ 850 nm) which is a strong absorption 
feature in the spectra. To compute the velocity dispersion pro¬ 
file we divided the pointing into six independent angular bins. 
In each bin, all spectra of all exposures where combined with 
a sigma clipping algorithm to remove any remaining cosmic 
rays. Velocity and velocity-dispersion profiles were computed 
using the penalized pixel-fittin g (pPXE) program developed by 
ICappellari & Emselleml (l2004l) . 

The SINEONI observations were carried out between 2008 
April and June (ESO proposal ID: 381 .D-0329(A); PI: Lanzoni), 
under an average seeing of ~ 0.8" (EWHM). The instrument was 
set to the 100 mas plate scale (pixel size: 0.05", 3.2" x3.2" field- 
of-view) using the K-band grating (1.95 - 2.45/im, R - 4000). 
The adaptive optics were performed using an R=12 magni¬ 
tude star located ~ 9" away from the cluster center. This re¬ 
sults in a Strehl ratio (the amount of light contained in the 
diffraction-limited core of the PSE, with respect to the total flux) 
of ~ 30%. Velocities were extracted individually for each star 
by matching the SINEONI pointing with high-resolution HST 
observations obtaining an accuracy of better than 0.2 spaxel. 
The spectra for the individual stars were extracted from the 
centroid positions of each star. The radial velocities were mea¬ 
sured fro m the CO band-heads using a Eourier cross-correlation 
method (iTonrv & DavisI Il979l) as implemented in the fxcor 
IRAE task. The velocity dispersion for the central region was 
then determined using these 52 individual veloc ities and us¬ 
ing the Maximum Likelihood method described in IWalker et al.l 
(1200^ . 

Eigure 11 in ILanzoni et al.l (1201 3l) clearly shows the large 
discrepancy between the two results from the different data 
sets. In their study they find a central velocity dispersion of 
CTp - 13.2 +1.3 kms“*, which is 40% lower than the result 
of the ARGUS observations (23.3 + 3.1 kms *). Both mea¬ 
surements, however seem to agre e wit h the outer kinematics ob¬ 
tained from ILanzoni et al.l (1201 3l) and iLapenna et al.l (1201 4l) us¬ 
ing GIRAEEE/ELAMES data. 

A concern is tha t the individual velocities as measured in 
ILanzoni et al.l (1201 3h may be biased towards the cluster mean, 
resulting in a lower velocity dispersion. We address this by mea¬ 
suring the velocity dispersion in magnitude bins from their pub- 
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SINFONI OBSERVATIONS SINFONI SIMULATIONS 



Fig. 2. Observed (left) and simulated (right) 3.2" x 3.2" field- 
of-view of the SINFONI IFU. The Strehl ratio of the simulated 
data cube is chosen to be the same as the observed one, 30%. 
Red crosses in the left panel mark the stars that were detected 
in the SINFONI observations, while green circles mark the stars 
that were not used in the final analysis of iLanzoni et al.l (1201 3h 
due to low signal-to-noise. 

ARGUS OBSERVATIONS ARGUS SIMULATIONS 


Fig. 3. Observed (left) and simulated (right) 7.2" x 8.1" field-of- 
view of the combined ARGUS IFU. The seeing of the simulated 
data cube is chosen to fit the observed data cube, 1.10". 


lished data. There should be no difference in the velocity dis¬ 
persion in these magnitude bins since the covered mass range 
is insignificant. If there is contamination from unresolved clus¬ 
ter light and light from the wings of nearby stars, we expect the 
effect to be larger for fainter stars resulting in the a lower dis¬ 
persion value. Figure[T]plots the velocity dispersion measured in 
three magnitude bins from their published data. It is clear that 
there is a trend to lower dispersion with magnitude signifying 
a potential bias. It is likely that even the brightest magnitude 
bin is also biased. By running a Spearman correlation test, us¬ 
ing the IDL routine r.correlate, on the absolute values of the 
individual velocities vs. magnitude we find a Spearman coeffi¬ 
cient ofR — -0.29 + 0.04 which indicates a slight anticorrelation 
with a two sided sign ificance of 2.1cr. Thus, it appears that the 
ILanzoni et al.l (1201 3ll central velocity dispersion is biased low. 
This trend raises a concern but is no proof for a bias. For this 
reason we set out to test for a possible bias in both datasets in 
detail. Contaminati on from neighboring and unresolved stars is 
studied in detail in iDubath et al.l (119941) . We explore a similar 
analysis to understand this bias, specific to NGC 6388. 


3. Simulated IFU observations 

We created a tool to simulate IFU observations for arbitrary 
weather conditions, instrument setups and objects. A detailed de¬ 
scription and verification of the code can be found in the upcom¬ 
ing paper Liitzgendorf et al. (2015, in prep.). Here we describe 
the concept of the code and its specific application to NGC 6388 
and the two observing techniques. 

Before creating the IFU dataset, a suitable star catalog is 
needed. To stay as close to NGC 6388 as possible we use the 
informatio n on the brightest s tars fr om a photometric catalog 
obtained in Liitzgendorf et al.l (1201 ll) from an image taken with 
the Advanced Camera for Surveys (ACS) on the Hubble Space 
Telescope (HST). We use an N-body realization to add the in¬ 
fluence of faint, unresolved stars to the catalog of NGC 6388. 
The model is set up by fitting an isochrone to the photometric 
catalog of NGC 6388. The best agreement is found with an age 
of 12 Gyr and a metallicity of z = 0.002 (Bressan et al., 2012; 
Chen et al., 2014). The isochrone is used to create an evolved 
mass function and to obtai n the stellar luminosity of individ¬ 
ual stars. Using an isotropic iKingl (1 1 9661) model, the distribution 
function with a central potential of Wo = 7 (values taken from 
McLaughlin et al., 2012) is integrated to retrieve the density pro¬ 
file p{r), the enclosed mass M{< r), and the potential <p{r). The 
stars are distribu ted in space using p{r) and a half-mass radius 
of r;, = 1.5 pc (iHarri ^ Il996h . As a final step, the distribution 
function and the potential are used to assign velocities to the 
individual stars in a way th at the cluster fulfills virial equilib¬ 
rium (see iHilker et aULOOTl) . We adjust the number of stars in 
the simulated cluster such that we obtain the same number of 
bright stars in the central pointing as measured in the observa¬ 
tions, to make sure that we do not overestimate the number of the 
fainter stars. The final cluster ma ss is 6.8 x 10^Mp, a bit lower 
than stated in the literature (~ l O^. lMcLaughlin & van der Marell 
12005^ lEutzgendorf et al.ll201ll) . To add the observed stars to the 
catalog, we interpolate their positions to the isochrone to obtain 
temperatures, surface gravities and K-magnitudes. We also re¬ 
move all stars from the catalog that are brighter than the faintest 
star from the SINFONI sample. In each step we assign new ve¬ 
locities to the stars by drawing randomly from a velocity distri¬ 
bution centred around 0 km s ' with a given velocity dispersion 
that we vary. 

In order to construct a realistic data cube from the simulated 
data we use a set of syn thetic spectra from the high-resolution 
synthetic stellar library (ICoelho et al.ll2005h obtained from the 
VIZIER arc hivj] for the ARGU S observations and synthetic 
spectra from iHusser et al.l (l20I3l) for the infrared observations 
with SINFONI. These libraries include spectra that cover the 
Calcium Triplet region that is used for the ARGUS observa¬ 
tions and the CO band heads to simulate SINFONI spectra. The 
spectra are first convolved with the spectral resolution of the re¬ 
spective instrument (ARGUS: R = 10400, SINFONI: R=4000) 
and re-sampled to the wavelength range to match a number of 
spectral elements of a typical observed spectrum from each in¬ 
strument. With this grid of parameters each star in the cluster is 
assigned a spectrum with a spectral shift accor ding to its radial 
velocity. As mentioned in ILanzoni et al.l (120131) . the unresolved 
background is in general featureless because of the weakening 
of the CO band heads with lower temperatures. This effect is 
also visible in our data, since the stars get assigned a spectra ac¬ 
cording to their temperature and surface gravity. However, for 


' Available at: http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=VI/120 
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completeness, we decided to keep the unresolved background 
even if it will not have a large effect on the measurements. 

The IFU construction is done by a C program that efficiently 
parallelizes the tasks of taking each star, finding a spectral match 
for it, taking its position and flux to construct a Moffat point 
spread function (PSF) 


f{x,y,a,f3) = {J3- 


1 + 



( 1 ) 


with p - 2.5 and a - FWHMI{2 V2F^ - 1), on the IFU grid 
and combines the spectra of all stars in each spaxel weighted by 
their luminosity. 

In the case of AO observations we use the combination of 
two Moffat functions, one with an intrinsic FWHM of 0.15" 
and one with the native seeing (i.e. the observed seeing, 0.8") 
weighted by the measured Strehl ratio. The construction of the 
final PSF for AO observations therefore follows the equation: 


f(x,y;ai,ao,P) ^ S * f(x,y;ai,P) + (I - S) * f(x,y;ao,P), (2) 

where S is the Strehl ratio and a, and Og are derived from 
the inner (intrinsic) and outer (native) FWHM, respectively. 
Furthermore, the routine applies noise to each spaxel. We request 
a minimum S/N of 10 for the faintest spaxel and apply noise to 
the remaining spaxels accordingly. FigureOcompares the central 
pointing of the SINFONI observations with the simulated data 
cube. The visual similarity of the two data sets gives confidence 
in the performance of the IFU simulation routine. The value for 
th e intrinsic FWHM is in good agreement with the values found 
in iRusli et al.l (l2013h . In addition, we select the 11 brightest iso¬ 
lated stars from the SINFONI cube and analyse their PSF. For 
each star we extract a ID PSF profile and fit a two-component 
Moffat function. The free parameters are the FWHM of the in¬ 
ner component (FWHM,), the FWHM of the outer component 
(FWHM„), the Strehl ratio (i.e. the ratio of these two compo¬ 
nents), and the scale. The fits show that the intrinsic FWHM does 
vary between 0.13" and 0.15". The Strehl ratios range from 0.1 
to 0.5, with the very brightest stars having the largest values. In 
general this shows that the parameters that we are using are a 
good representation of the observed data. Furthermore, we di¬ 
rectly compare the ID profiles of our simulated SINFONI data 
cube and the observed one and find very good agreement with 
the chosen values. 

The simulated PSFs assume two simplifications: Spherical 
symmetry and no variability across the field of view. The ellip- 
ticity that is observed in the measured datacube is difficult to re¬ 
produce in the simulations and would introduce even more free 
parameters. We concluded, however, that the effect of the elliptic 
shape of the PSF will be small to the overall analysis. The ellip- 
ticity will cause some of the stars (that are aligned with the minor 
axis of the ellipse) to be less contaminated by bright stars while 
other stars (aligned with the major axis) will be more contami¬ 
nated. This would have an effect if we would make a star-by-star 
analysis. However, since we are computing the velocity disper¬ 
sion of the system, these effects will cancel each other out. The 
same argument holds for possible variability of the PSF across 
the field of view. The ID fits of the PSF show that our models 
are compatible with the observations and that the fitted PSF pa¬ 
rameters scatter around the parameters that we use for the final 
simulations. 


4. Extraction of kinematics 

The extraction of the kinematics is performed in a similar way 
as with the observed data. For SINFONI we use the star position 
from the catalog to extract the spectrum of the central spaxel of 
each star from the IFU. We note that by skipping the process 
of cross correlating with an independent catalog we neglect the 
errors that might arise from mismatches with the catalogs but 
we consider them as marginal since the IFU data cube is con¬ 
structed using the same catalog.From the extracted spectra we 
then measure the velocities of each star using fxcor from the 
IRAF package and a list of templates with the same stellar pa¬ 
rameters as the stars. We make the assumption that the internal 
parameters of the stars are perfectly known and that we are using 
the perfect template for each star. The final velocity dispersion 
is then com puted by using the Max imum Likelihood method in¬ 
troduced bv iPrvor & Mevlanl(ll993h to account for the individual 
uncertainties on the velocity measurements. We repeat this mea¬ 
surement for IFU simulations with different Strehl ratios ranging 
from 10% to 100%. The results are discussed in the following 
section. 

For the simulated ARGUS observations we adopt the po¬ 
sition of the center and the radial bins from lLiitzgendorf et al.l 
(1201 ih . We then combine all spectra in each bin by applying an 
iterative sigma clipping scheme. The velocity dispersion is de¬ 
rived by performing fits with pPXF on the combined spectrum 
but this time not measuring the centroid (i.e. the velocity) but 
directly the velocity dispersion from the broadening of the lines. 
For both SINFONI and ARGUS simulations we repeat the pro¬ 
cedure described above for 100 iterations in which we assign 
new velocities to all the stars in each step. This is performed for 
four different input velocity dispersions (10, 15, 20, 25 kms^') 
and 10 different values for seeing and Strehl ratios. The final 
values and their uncertainties are taken from the mean and the 
standard deviation of all iterations. 


5. Results 

5.1. Velocity measurements 

Figure |4] shows the result of the two IFU simulations. For the 
SINFONI observations we show the obtained velocity disper¬ 
sion as a function of Strehl ratio and for different input velocity 
dispersions. For the ARGUS observations the same is shown but 
for different values of seeing. The Figure shows that despite large 
uncertainties, the ARGUS measurements closely resemble the 
input velocity dispersions (dashed lines of the same color). The 
SINFONI observations, however, are strongly biased towards 
lower velocity dispersions. This can be explained by individual 
and especially faint stars being contaminated by the wings of 
their neighbouring stars. This pushes measured velocities closer 
to the mean velocity of the cluster and therefore lowers the ve¬ 
locity dispersion. As shown in Figure |4] this effect is especially 
severe for low Strehl ratios but even with a perfect Strehl ratio 
of 100% the bias persists. We find that a measured veloc i ty dis- 
persion of about 13 km s ' (as measured in lLanzoni et al.l (1201 3h 
with a Strehl ratio of 30%) requires an input velocity dispersion 
of ~ 25 km s '. T his is in excellent agreeme nt with the measure¬ 
ments obtained in iLiitzgendorf et al.| (1201 ll) and would explain 
the discrepancy between the two measurements. 
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Fig. 4. IFU simulations for SINFONI and ARGUS (central bin) with different input velocity dispersions as a function of Strehl ratio 
and seeing. The input velocity dispersions are shown as dashed lines, the measured velocity dispersions as dots in the corresponding 
color. 
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Fig. 5. Same as Figure|4]while fixing the velocities of the three brightest stars to their measured velocities and different input velocity 
dispersions as a function of Strehl ratio and seeing. 


5.2. Discussion of the bias 

The bias towards lower velocity dispersion in the SINFONI 
dataset is concerning. To investigate further, we plot the mea¬ 
sured velocities versus the input velocities of the SINFONI stars 
for a set of 10 Monte Carlo runs. We also color code the stars by 
their magnitudes. Figure shows the measured velocities. It is 
obvious that the bias depends on the magnitude of each star and 
that for the very bright stars (dark red - red) the velocities are 
measured well, while for the faint stars (blue) the bias is always 
directed to a lower velocity. 

The reason for this bias can not be explained by the unre¬ 
solved background, since, as discussed in the previous sections, 
this is featureless. However, the influence of the neighboring 
stars and the contamination of the faint stars by the wings of 
the bright stars must not be underestimated. To quantify this 
effect we compute a contamination map for each spaxel in the 
SINFONI pointing. Figure [T] shows the maximum contribution 
to each spaxel by any star in the left panel. From this map we 
created the contamination map in the right panel where we ex¬ 
amined the spaxel of the measured stars in the SINFONI point- 
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Fig. 6. Measured velocities versus input velocities for 100 Monte 
Carlo runs on the SINFONI dataset. The stars are color coded 
according to their magnitudes. 
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Fig. 7. Contamination map for the SINFONI pointing. Left panel shows the maximum contribution of a single star to each spaxel. 
From this map the right panel is constructed that shows for each star that is measured in the SINFONI pointing the percentage of 
contamination by other stars. All stars above the dark blue regime are contaminated by more than 20%. 


ing and determined to which level these spaxel (that were used 
to derive the velocity) are contaminated. To give an example we 
selected two stars close to the brightest star (A) with different 
magnitudes and state their contamination. While star B is already 
moderately contaminated (30%) by the bright star and other stars 
around it, star C is highly contaminated by the wing of the very 
bright star A by a value of ~ 60%, i.e. it contributes less than 
half of the light to this spaxel. This shows the importance of the 
PSF wings even in AO observations and the contaminating ef¬ 
fect of the bright stars. For a Gaussian velocity distribution, the 
contamination from the neighboring stars will bias the velocities 
of the faint stars to the mean velocity. We note that we also re¬ 
run the simulations using a smaller value for the outer FWHM 
(FWHM„ = 0.5") in order to test the effect of the native PSF. 
We find the bias reduced by 10 - 15% but still present and ar¬ 
gue that a seeing of 0.5" is not representative of average VLT 
observations and much lower than the observed seeing of 0.8". 

Compared to the SINFONI data set, the simulated ARGUS 
data in the right panel of Figure |4] does not show such a bias. 
However, the simulations of the innermost bin result in a very 
large scatter and an uncertainty of more than 20 %. Similar re¬ 
sults have been already reported by Bianchini et al. (2015, sub¬ 
mitted). The highest values for cr also show a slight trend to¬ 
wards lower values, but still consistent with the input velocity 
dispersion, when taking the large uncertainties into account. This 
lower values most probably arise from the fact that the innermost 
bin is dominated by one or two single stars with narrow absorp¬ 
tion lines. Not shown in Figure|4]are the results for the five other 
bins of the ARGUS IFU that were also included in the simula¬ 
tions. The resulting uncertainties are smaller than for the central 
bin and are shown in Figure [8] 

In order to follow the argument bv iLanzoni et al.l (1201 3h that 
the brightest stars might bias the velocity dispersion towards 
higher values for this specific case of velocity distributions we 
run the same simulations as before but this time with fixing the 
velocities of the three brightest stars to the measured veloci¬ 
ties of vi = -23.2 kms^*, V 2 = 18.4 kms^*, (measured by 


iLanzoni et akl (l2013h . the two brightest stars in the sample) and 
V 3 = -46.6 k m s~* (measured from the ARGUS pointing and not 
present in the ILanzoni et al.l (1201 3l) dataset. The velocities of the 
rest of the stars are again drawn from a velocity distribution with 
different input velocity dispersion. Figure |5] shows the result of 
these runs. While the SINFONI dataset is only slightly affected 
(compared to Figure |4]i, the impact of the three bright stars hav¬ 
ing opposite velocities is visible in the ARGUS simulation. For 
low velocity dispersions the three stars drive the velocity dis¬ 
persion up to 20 km s ' but for high velocity dispersions there is 
barely a difference (while the bias in the SINFONI data is visible 
for all velocity dispersions). The reason for this directed bias lies 
in the fact that adding two very fast stars on top of a velocity dis¬ 
tribution with a low velocity dispersion for all simulation runs is 
not a correct representation of this distribution. Further, we want 
to stress that a velocity dispersion as low as 18 kms“* is still al¬ 
lowed by the uncertainties. As a last point it is worth mentioning 
that only the central point is heavily affected by this specific con¬ 
stellation of stars but not the bins further out. As we show in the 
next section, only the bins further out have an effect on the black 
hole mass in the Jeans modelling, due to the large errorbars on 
the central point. 


Thus, we conclude that even excluding the central region 
where the bright stars are, there still is a rise in the dispersion 
at larger radii where these bright stars have no influence. The 
observed rise in the dispersion is therefore not dependent on 
the bright stars. Second, one cannot remove the stars from the 
sample, unless there are reasons due to them not being in dy¬ 
namical equilibrium. Our simulations are designed to measure 
the uncertainties on the dispersion, including the brightness of 
the stars. We must include them when measuring the dispersion, 
and then we determine the uncertainty from the simulations. It is 
irrelevant what their velocity distribution is for the simulations. 
Furthermore, we believe that the possible rotation signature in 
the center of NGC 6388 is not due to three bright star s but t o 
a group of stars as already stated in iLfitzgendorf et ak (201 ll). 
The reason why this is not seen in the dataset of Lanzoni et al.l 
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Fig. 8. Jeans models for the velocity dispersion profile o f NGC 
6388 using the profile obtained in iLtitzgendorf et alJ (1201 ih 
(dots), the uncertainties on these measurements obtaine d in this 
work, the corrected value from iLanzoni^t^ J 2()T^ It rian- 
gle), and the outer po ints obtained by Lanzoni et al.l (1201^ and 
iLapenna et al.l (1201 4 squares). 


(1201 3h is assumed to be due to the bias of the fainter stars to¬ 
wards the mean velocity which leads to a weakening of the rota¬ 
tion amplitude. While any possible rotation is interesting for the 
cluster kinematics, it must be included for a dynamical analy¬ 
sis for the enclosed mass. The second moment we measure from 
the integrated light naturally includes both the ordered and hot 
components of the stellar kinematics. 

5.3. New dynamical models 

As mentioned in the previous paragraph, we simulate the re¬ 
maining five bins of the ARGUS observations in the same 
way as described abov e and find that the uncertainties in 
iLutzgendorf et al.|(l201 ih are underestimated (see Figure[8]l. The 
difference might arise from th e more advanced IFU mo delling 
compared to the method used in lLiitzgendorf et ^ (1201 Ih where, 
for example, only one spectral type of star w as used. We there¬ 
fore re run our Jeans models as described in ILtitzgendorf et al.l 
(1201 111 , using the same, but slightly smoothed surface brightness 
profile in addition with a M/L profile obtained from N-body sim¬ 
ulations shown in Figure|9] The M/L profile is obtained by fitting 
a grid of N-body models to the photometric and kinematic data 
of NGC 6388 and computing the M/L for the best-fit model. 
The models are computed with 10% stellar-m ass black holes 
and neutron stars retention fraction and an initial IKt^ HidMl 
mass fun ction. More detai l s abou t the N-body models are pre¬ 
sented in lMcNamara et all (l2012l) . 

For completeness we also add the m ost re ce ntly published 
outer datatpoints from iLanzoni et al.l (l2013h : iLa^ennaet^ 

( 20JA , squares) and the corrected datapoint from ILanzoni et al* 
(1201 3l) (triangle, correction is indicated with the dashed arrow, 
not included in the fit). Figure [8] shows Jeans models with differ¬ 
ent black-holes masses, the projected M/L profile and its MGE 
parametrization, as well as the values for each model are 
shown in Fi gure E We no t e that our Jeans models differ from 
the ones of ILanzoni et al.l (1201 3l) by the M/L profile and re¬ 
sult in a higher model central velocity dispersi on even with- 
out a black hole. Furthermore, models used by ILanzoni et al.l 


(12013h are the re sults of parametric multi-mass iKingl (Il966h 
and WilsonI (Il975h models while our models only take the mea¬ 
sured and smoothed surface brightness profile and a M/L pro¬ 
file as input.Furhtermore, it should be stated that the inner¬ 
most datapoint does not influence the black-hole measurement 
due to its large uncertainty. However, the bins further out still 
show the clear rise and cannot be the result of contamination 
by the innermost bright stars. The models result in a black- 
hole mass of M, - (2.8 + 0.4) x 1Q‘^Mq and a total MjL - 
(1.6 + 0.1)Mo/Lo. We note that the uncertainties on the black- 
hole mass might be underestimated in this case since we do not 
perform Monte Carlo simula t ions o n the surface brightness pro¬ 
file as in ILtitzgendorf et al.l (1201 ll) . Both the black-hole mass 
and t he M/L agree well within their uncertainties with the results 
from lLiitzgendorf et al.l(l201 ll) (M. = (1.7+0.9)xlO"^Mo, M/L = 
(L6±0.3)Mo/Lo). 


We note that we are using simplified models to obtain 
the black-hole mass. By using isotropic Jeans models we are 
neglecting the fact that a globul ar cluster most likely sh ows 
anisotropy in its outer regions. In lLiitzgendorf et aP (1201 ll) we 
demonstrated that in the central regions of a globular cluster 
anisotropy quickly disappears due to relaxation processes. By 
fitting the entire cluster profile, ho wever, we include r egions 
that can be affected by anisotropy (IZocchi et alJ 1201^ . This 


More sophisticated models such as 
20151) or Schwarzschild Models (e. 


f, (e.g. Zocchi et alJ 

2012 

Ivan den Bosch et al.l 

2006 


Jardel & Gebhardi2013l) would therefore be the better represen¬ 


tation of such a complex system, but are outside of the scope of 
this paper. 


6. Summary 

We perform IFTJ simulations of observations with two different 
IFUs for the particular case of the globular cluster NGC 6388. 
Using a combination of synthetic spectra, N-body realizations 
and observational data we reproduce observations of the central 
region of NGC 6388 as close to reality as possible. The kine¬ 
matics are extracted in similar ways as done for the observed 
data. That allows a direct comparison of the biases coming with 
the technique of measuring individual velocities and line broad¬ 
ening of integrated light. The results are evaluated for different 
observing conditions and input velocity dispersions. 

We show that both measurements feature large uncertainties 
but only the SINFONI datasets are clearly biased towards lower 
velocity dispersions depending on the Strehl ratio of the observa¬ 
tions. The bias towards lower velocities can be explained by the 
contamination of neighbouring stars which shifts the velocities 
of the individual stars closer to the mean velocity. For a mea¬ 
sured velocity dispersion of about 13 kms“', an input velocity 
dispersion of about 25 kms"' is required. Thi s is in good agree¬ 
ment with our measurements presented in ILtitzgendorf et al.l 
(1201 ih . In the integrated light method on the other hand, stars 
and background light are both valuable signal that is used to 
derive the velocity dispersion. However, we noticed a bias to¬ 
wards higher velocity dispersions using integrated light when 
fixing the extreme velocities of the three brightest stars regard¬ 
less of the input velocity dispersion. The bias is included in 
the large uncertainties of the innermost datapoint and disap¬ 
pears for high velocity dispersions. We run Jeans models on 
th e whole kinematic da taset including the o uter p oints derived 
in ILanzoni et al.l (1201 3h and iLanenna et al.l (l2014l) . and includ¬ 
ing the newly derived uncertainties for our points measured 
with the ARGUS IFU. The fits result in a black-hole mass of 
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Fig. 9. Input M/L profile obtained from N-body simulations (left) and;^f^ for the different Jeans models from Figure [8] (right) with a 
best-fit black-hole mass of M, - (2.4 + 0.4) x 10 '*Mq. While the rise in the M/L in the central region is due to stellar remnants, the 
rise at larger radii originates from low-mass stars. 


M. = (2.8+0.4)xl0'*Mo and M/L ratio M/L = (L6+0.1 )Mq/Lo 
which is in good agreement with the previous results on this data 
set. 

This analysis is specific for NGC 6388 where there are many 
observations and photometric data to compare with. However, 
biases for measurements with individual velocities are most 
likely present in all of the data sets that have been obtained so 
far. NGC 6388 belongs to the densest globular clusters in our 
Milky Way. Therefore it is understandable that the blend effects 
are severe in this environment. However, even with sparse clus¬ 
ters one has to be careful when taking individual velocities and 
tests for their reliability are always needed. For core-collapsed 
clusters with densities even higher than NGC 3688, it is not rec¬ 
ommended to use this technique at all. One solution to the is- 
su e might arise from a n ew deconvolution technique developed 
bv iKamann et al.l (l2013l) . The precise extraction method reduces 
the overblend effects of neighboring and background stars to a 
minimum and allows to analyze data with much lower spectral 
resolution than the integrated method. 

For the case of NGC 6388 and other nearby clusters there 
are other methods for measuring the velocity dispersion. Proper 
motions for a large sample of Galactic globular clusters have 
been obtaine d in the frame of the A CS survey of Galactic glob¬ 
ular clusters (ISaraiedini et al.ll2007h . Usin g proper motions to- 
gether with discrete Jeans modelling (e.g. IWatkins et al]l2013h 
brings a vast advantage to globular cluster modelling. By in¬ 
cluding M/L profiles and anisotropy, these models do not just 
apply the physics needed to the models but al s o do n ot need to 
make assumptions for binning. IWatkins et alJ (12015h published 
the kinematic profiles of a large set of globular clusters includ¬ 
ing NGC 6388. Unfortunately their profile does not extend fur¬ 
ther inside than 1 arcsec and therefore does not cover the critical 
region that is discussed in this paper. If there will be enough us¬ 
able data points on the central parts of NGC 6388 measured in 
the future, proper motions might help to solve the discrepancies 
in radial velocity measurements. 

A detailed and more general study on the influence of IFU 
measurements for semi-resolved objects such as globular clus¬ 
ters as a function of cluster properties and observing conditions 
is crucial to produce reliable results. In terms of non-AO sup¬ 
ported observations similar to the ARGUS data set and the de¬ 


convolution extraction technique by iKamann et al.l (12013h we 
have performed detailed studies on the reliability of velocity dis¬ 
persion measurements with IFUs and their effect on the determi¬ 
nation of the total mass and the black hole mass in the center, 
if present. The results of this study will soon be published in 
Liitzgendorf et al. (2015, in prep.). This provides a valuable tool 
for future and present observations and a key to calibrate mea¬ 
surements performed with IFU instruments. 
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